Numerical Analysis of Optimized Coherent Control Pulses 
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Numerically we simulate the effect of optimized coherent control pulses with a finite duration on a 
qubit in a bath of spins. The pulses of finite duration are compared with ideal instantaneous pulses. 
In particular, we show that properly designed short pulses can approximate ideal instantaneous 
pulses up to a certain order in the shortness of the pulse. We provide examples of such pulses, 
quantify the discrepancy from the ideal case and compare their effect for various ranges of the 
coupling constants. 
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I. INTRODUCTION 

The coherent control of quantum systems continues to 
be a topic of great interest. The possibility of maintain- 
ing a spin in a coherent state is of extreme importance 
in fields of application like nuclear magnetic resonance 
(NMR) or the manipulation of quantum dots. In particu- 
lar for quantum information processing, a long coherence 
time of the qubit is an indispensable prerequisite for its 
realization. 

A quantum bit (henceforth: qubit) is a two-level sys- 
tem which is conveniently regarded as a spin S = 1/2. 
Operations on qubits to change or to correct their state 
are performed through quantum gates. Their effect on 
the density matrix of the qubit can be described as a rota- 
tion in the Bloch sphere. Experimentally, they often can 
be obtained by the application of electromagnetic pulses. 
A 1-qubit gate is generally a single rotation about a given 
axis a in spin space. The angle of rotation classifies the 
type of the pulse. For instance a tt pulses rotates the spin 
by 180°. These pulses find a wide range of applications 
in dynamical decoupling [J, @, H, 0, 1, S 0| and in NMR 
d, Q where also tt/2 pulses are crucial. In quantum in- 
formation processing the 7r/2 pulse in combination with 
a 7r pulse realizes the important Hadamard gate. 

The idea of dynamical decoupling (DD) 0, H i, 0, i, 
H, 0] has been developed from the spin echo technique 
in NMR [H [H [3. DD aims at decoupling the qubit 
from the environment by means of the application of ap- 
propriate pulse sequences. From a theoretical point of 
view the topic has been widely studied and many dif- 
ferent sequences of pulses have been proposed. Among 
these we recall the series of periodic equidistant 7r pulses, 
called bang-bang control (BB) |2j], the periodically it- 
erated 2-pulse sequence accor ding to Carr/Purcell and 
Meiboom/Gill (CPMG) [1, HH, E3 ; the concatenated se- 
quence (CDD) proposed by Khodjasteh and Lidar [H,[l3| 
as well as the fully optimized sequence (UDD) derived by 
one of the authors @, 0j| . 

Experimentally, the spin echo and the CPMG sequence 
are standard in NMR j8[. To our knowledge, other se- 
quences have not yet been tested. In realizations of qubits 



on the basis of semiconductor technology so far only the 
spin echo technique has been implemented 0, [lT], EH • 
But computations for quantum dot systems show that 
more elaborate pulse sequences are very likely to be use- 
ful in suppressing decoherence, see for instance lol Rj [Til . 

Most theoretical examples (for exceptions see Ref. 
[HI) reported so far have the limitation that the pulses 
are assumed to be ideal. This means that the pulses are 
considered to be instantaneous and infinitely strong in 
the sense of a 6 peak. In this case, one is allowed to 
ignore the effect of the bath, inducing the decoherence, 
during the action of the pulse because the coupling to 
the bath is negligible relative to the amplitude of the S 
pulse. Hence the rotation due to the pulse can be viewed 
to be completely separate from the free evolution of the 
system, qubit and bath, without pulse. 

If the pulse has a finite duration (r p ) so that its time 
of application is comparable with the characteristic time 
scales of the bath, the separation between evolution due 
to the pulse and evolution of the undriven system is not 
valid anymore. If we suppose that the duration t p is still 
small an expansion in t p about the limit of a 5 pulse 
is appropriate. The proposed scenario [2(| establishes 
an equivalence, up to corrections expanded in a series in 
r p , between the real pulse and an ideal S pulse at some 
intermediate instant t s with < r s < t p , see Figs. [TJ [21 
and [31 Before and after the ersatz pulse at r s , the free 
evolution of the system, qubit and bath, without pulse 
takes place. 

The corrections expanded in powers of t p depend also 
on the shape of the pulse; so one can aim at making them 
vanish or at minimizing them by shaping the pulses skill- 
fully. This is the route that we established previously [2(| 
analytically by the expansion in r p . In the present work, 
we demonstrate numerically that the higher order correc- 
tions neglected in the analytical calculations are indeed 
negligible. Thereby, we have shown not only the valid- 
ity of the previous analytic calculation but we have also 
demonstrated that the real performance of the proposed 
pulses is advantageous. 

We draw the readers' attention to the fact that shaped 
pulses have been introduced in NMR previously, see for 
instance the Refs. M, II [H H M, M and Ref. © 
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Figure 1: (Color online) Examples of it pulses implemented 
in the simulations. The ideal pulse is given by a S peak op- 
erating at the instant r s . SGLPi is the standard pulse of 
constant amplitude without optimization of the pulse shape. 
For details see Table U 
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Figure 2: (Color online) Examples of 7r/2 pulses implemented 
in the simulations. The ideal pulse is given by a 8 peak op- 
erating at the instant r s . SGLPi2 is the standard pulse of 
constant amplitude without optimization of the pulse shape. 
For details see Table ITU 



for an overview in the field of quantum information. But 
the goals of these investigations were different from ours 
even though it turned out that for tt pulses certain shapes 
with t s — Tp/2 happen to coincide [20| . 

The paper is organized as follows. In Sect. [IT] we briefly 
recall the analytical arguments for the expansion in pow- 
ers of t p ; especially the expressions for the first and sec- 
ond order corrections are given. Then we introduce a 
quantity to measure the deviation of the real pulse from 
the ideal pulse and compute this deviation analytically. 
Moreover, we relate the two parameters of the model 
to the experimental situation in various realizations of 
qubits. In Sect. IHII the spin Hamiltonian is introduced 



Figure 3: (Color online) Examples of n/2 pulses for which also 
some of the second order corrections vanish, namely 7721 = 
and 7722 = (see main text). For details see Table ITT1 



which serves as our system of a qubit coupled to a deco- 
herence bath. For this model we compute the deviation 
between real and ideal pulse analytically and numerically. 
The experimentally relevant ranges of parameters are es- 
timated. The numerical results are discussed in Sect. IIVI 
for 7T and for ir/2 pulses. Finally, in Sect. [V] we draw our 
conclusions. 



II. THEORETICAL PREDICTIONS 

A. First and second order corrections 

In order to disentangle the actual pulse and the free 
evolution of the system we proceed as follows. The total 
unitary time evolution during the real pulse is split into 
the time evolution of the system alone and of the pulse 
alone which is taken to occur at t s within the interval 
[0,t p ], see Ref. 20. The time evolution of the system 
alone is taken to occur before and after the evolution 
due to the pulse. The evolution due to the pulse is mul- 
tiplied additionally by corrections coming from the non- 
commutation of the Hamiltonians of the pulse and of the 
system. They can be expanded in a series in r p . It is 
important to stress that this technique does not aim at 
eliminating the coupling between the qubit and the bath 
completely, but only at separating the effect of the pulse 
from that of bath. The coupling between the qubit and 
the bath remains active during the free evolution of the 
system. 

To be explicit, we consider the following general Hamil- 
tonian 

H tot = H + H {t), (1) 
where the Hamiltonian H of the qubit coupled to the 
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bath is 



H = H b + XAcr z , 



(2) 



where H b is a completely general bath and A a completely 
general coupling operator acting on the bath. The Pauli 
matrices represent operators acting on the qubit. The 
internal energy scale of H b shall be denoted by ui b while 
A is the coupling constant between qubit and bath. 

Note that we assume only a coupling along the z di- 
rection. Hence the model contains only dephasing, i.e., 
a finite T 2 . No spin flips are possible so that T\ = 00. 
Though this represents a restriction it is well justified for 
large magnetic fields along z so that all other couplings 
average out in the rotating-frame approximation. 

The Hamiltonian of the pulse is denoted by Hq 



H (t) = v(t)c. 



representing a rotation around the y axis. The pulse 
shape is given by the function v(t). Note that Hq and H 
do not commute implying that the unitary time evolution 
U(t p , 0) during the application of a pulse is a non-trivial 
quantity. 

Splitting the time evolution U(r p , 0) into the time evo- 
lutions during two intervals, U(t p ,t s ) and U(t s ,0), and 
formally solving the Schrodinger equation for each of 
them with a suitable ansatz we eventually obtain (for 
details see Refs. uO) 



U p (t p ,0) = T{e-^o Tp ff*o t (t)dt| 

= e -i(T p -T B )H e -i*y ft <W UF ( Tp , 0) 



e ~ ia y Jo" v(t)dt e -iT 3 H _ 



(4) 



where Uf{t p ,0) represents the correction term. Without 
any correction, i.e., for Uf(t p ,0) = 1, the two exponen- 
tials of the pulse can be combined in the middle of the 
right hand side of Eq. (|4]) so that the unitary operator of 
the ideal pulse occurs 

U p (T p , 0) = e -^v-r s )H e - lay JJ» v(t)dt e -ir 3 H^ ( 5 ) 

The correction is expanded in a series in powers of t p 
U f (t p , 0) = exp(-z(77« + V W + ...)) (6) 
where r/ J ) is the term of order t^ . We obtained [2(| 

= (Vii<?x +m2cr z )XA (7a) 
i {miVx + V220- z ) X[H b ,A] + r\ 2Z a y X 2 A 2 . (7b) 



,(2) _ 



Note that [Hi,, A] is of the order of ui b so that the corre- 
sponding term is indeed of order Xu> b T 2 , thus of second 
order in t 2 . This becomes manifest in the explicit inte- 
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-0.00653 
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0.18087 
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Table I: Overview of the n pulses implemented in the simula- 
tions. UPi and SGLPi are symmetric pulses (r s = t p /2). The 
switching instants r, and the amplitudes are given in units of 
t p and 1/t p , respectively. The column r/ 2 ) refers from top to 
bottom to the coefficients 7721 , 7722, 7723 in units of r p , see Eq. 



gral equations for the coefficients rjij 
Vu = ( t p - T s) sin ip Tp + t s sin i/j - / sin rp t dt 
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'(t p -t s ) 



am(ipt — ipo) dt 
sin(i/v p - ipt) dt 



sin(?/) tl - ■0t 2 )sgn(*i - t 2 )dt 1 dt 2 , (8e) 



where ipt — 2 J v(t')dt', At = t — T s , and 8 = t/j Tp — ipo is 
the area under the amplitude of the pulse. The angle 8 
represents the total angle of rotation of the qubit 's spin 
under the action of the pulse. 

The function v(t) and the instant t s are the free vari- 
ables which can be fine-tuned to ideally make the coef- 
ficients rjij vanish or at least to minimize their moduli. 
In Fig. Q] examples of piecewise constant pulses for 8 = it 
are reported. The pulse SGLPi is the standard pulse of 
constant amplitude which has finite first and second or- 
der corrections. The pulses UPi and ASYPi are chosen 
such that their first order correction jjW vanishes. Their 
second order correction rj^ does not vanish. We have 



proven previously that rf 2 ^ cannot be made vanish for a 
7r pulse. [2(|. For quantitative details, see Tab.Q] 





amplitude(s) 


n 




SGLPi2 


1/2 


tt/4 






UPi2 


1/2 


±1.65765 


0.13155 
0.86845 


-0.01305 


0.05151 


ASYPi2 


0.23128 


±1.39116 


0.78220 


-0.01279 
-0.05691 
0.88990 


S2NDPi2 


1/2 


±2.31993 


0.05848 
0.22384 
0.77616 
0.94152 





±0.01335 


A2NDPi2 


0.61218 


±2.09429 


0.08361 
0.29828 
0.90217 





±0.01659 



Table II: Overview of the tt/2 pulses implemented in the sim- 
ulation. UPi2, SGLPi2, and S2NDPi2 are symmetric pulses 
(t s = t p /2). The switching instants n and the amplitudes 
are given in units of t p and 1/t p , respectively. The column 
j/ 2 ' refers from top to bottom to the coefficients 7721 , V22 , 1)23 
in units of Tp, see Eq. ([7]) 

In analogy, Fig. [5] depicts examples of piecewise con- 
stant pulses for 9 = tt/2. The pulse SGLPi2 is the 
standard pulse of constant amplitude which has finite 
first and second order corrections. The pulses UPi2 and 
ASYPi2 are chosen such that their first order correction 
t/ 1 ) vanishes. Their second order correction ?/ 2 ) does not 
vanish. For the quantitative details, we refer the reader 
to Tab. OH 

The pulses S2ND2 and A2NDPi2 are plotted in Fig. 
121 They are chosen such that their first order correction 
t/ 1 ) and the second order coefficients 7721 and 7722 vanish. 
We were not able to find a solution which has addition- 
ally r/23 — 0. But we have not succeeded in proving the 
impossibility of finding such a solution either. For the 
quantitative details, we refer the reader to Tab. HT1 
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Figure 4: (Color online) Sketch of the spin chain representing 
qubit and spin bath. The qubit 1 is coupled to the spin 2 of 
the chain. The coupling between qubit and the spin bath is 
given by J while the internal exchange coupling within the 
chain is a J. 

the ideal time evolution (Up) from the real one (Up). To 
capture this distance by a single number we define the 
norm 



d := ymax{Eigenvalues(AtA)}. (9) 
For a pulse of angle 6, the ideal pulse reads 

Up = e -i(r p -T 3 )B e -i^y e -ir 3 H ^ 

while the real pulse is given by 

U*p = e-^- T - 3ff e- i ^ <7 "C/ i r(T P! 0)e i ^ 1 ^e- iTsK , (11) 
such that 

A = e - i(T >>- T ° )H e- i± ^' T »(I - U F (T p ,0))e l ^^e~ lT " H . 
This implies 

At A = E/t(l - U f (t p , 0))(1 - U f (t p ,0))U, (13) 

where U — e~ l ^ L ' Ty e~ 1TsH is a unitary operator which 
leaves the eigenvalues of the product (1 — U f (t p , 0))(1 — 
Uf(t p ,0)) unaffected. 

Next, we expand Uf in r p . If the leading order is rf 1 ' 
we have Uf ~ 1 — vrf^ + C( r p) whence 



d = \ / max {Eigenvalues(?7( 1 )?7( 1 ))} ± 0(t|). (14) 

If the leading order is -q^ we have Uf ~ 1 — irf^ +0(Tp ) 
whence 

d = y'max {Eigenvalues^ 2 )?/ 2 ))} ± 0(t^). (15) 

We deduce that in the case of finite first order t/ 1 ) ^ 
one has d = 0(t p ) while for vanishing first order, but 
finite second order rf 2 ' 7^ one has d = 0(ri). 



B. Measure of deviation 

The above results represent the analytical finding that 
we intend to check numerically. In order to do so we 
need a measure of how well the real pulse approximates 
the ideal instantaneous one. We define the operator dif- 
ference A := Up — Up which quantifies the distance of 



III. THE SPIN CHAIN AS DECOHERENCE 
BATH 

A. The Model 

The formulae (|14ll5p for d hold for any Hamiltonian 
that can be expressed in the form Next, we specify 
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the model we investigate numerically. It is a spin chain 
of N spins where the first spin (a) represents the qubit, 
see Fig. |4] The Hamiltonian considered is given by 

N 

H s = Ju z uf +ajJ2 ■ ■ (16) 

i=2 

Obviously, this Hamiltonian is an example for the most 
general dephasing Hamiltonian ([2]) . In (p~6|) the bath is a 
bath of spins and the coupling between bath and qubit 
is quantified by J; hence we have A = J. The internal 
energy scale of the bath uif, equals a J in (Til)]) . 

In order to apply our general results (|14I15[) we have to 
compute r/ 2 ) for the specific case of the spin Hamiltonian 
(IT5|) . The bath operator A in (J2J) consists only of the z 
component of the second spin. Hence one has A 2 = 1. 

The other term in (|7bp comprises [Hb, A]. For (fl6j) this 
commutator contains only the second and the third spin. 
Hence we anticipate that the numerical results will not 
show any significant size dependence in the regime where 
the expansion in t p is valid, i.e., for low values of A and 
Wb which translates to low values of J. Explicitly we find 
for rjW 

7? (2) = -2J 2 a(r ]21 a x + r l22 a z )(a (2) x o t(3) ) 2 + 7 ]23 J%, 

(17) 

where ( ) z stands for the z component. Because only 
three spins occur, it is a basic exercise to determine for 
r/ 2 ** given by (I17|) the maximum eigenvalues of (tj^) 2 
yielding 

d = J 2 ^16 a 2 ( V 2 21 + V 2 2 ) + V 2 3 + 0(t5). (18) 

In this formula, the quadratic dependence of d as a func- 
tion of J has been put in evidence. The quadratic depen- 
dence on t p is less manifest, but it becomes obvious on 
inspecting the integrals in {8} from which r)2j — 0(t 2 ) 
ensues. 

Once t s and v(t) are known, the coefficients rj 2 i, 
7/22, and r]23 can be easily computed according to (jSJ. 
Thereby, we have an analytical prediction for the lead- 
ing order of d as function of J including the prefactor. 
For fixed value of </, Eq. (| 18[) as a function of a is char- 
acterized by a constant behaviour dominated by 7723 for 
a <C 1 and a linear behaviour in a for large values of the 
coupling constant. 

B. The Range of Parameters 

Although we are focusing here on the theoretical issues 
it is helpful to have an idea about the experimental range 
of parameters. In the sequel, we thus try to assess the 
relevant ranges. The numbers given represent only crude 
estimates since the precise values depend strongly on the 
particular experimental setup. Moreover, the relevant 
decoherence processes are not yet always known. 

First, we consider liquid NMR like crotonic acid or 
alanine [28| . The pulse lengths r p used are in the range of 



200/js. The maximum pulse amplitude B m for a it pulse 
is thus in the range of 10kHz. The couplings between the 
nuclear spins lie between 1 and about 70 Hz. A key ratio 
is J/B m , i.e., the relative dimensionless strength of the 
pulse. Here it takes values in the range of 10~ 4 and 10~ 2 . 
The other important parameter a is the dimensionless 
ratio uib/X between the internal energy scale u>b of the 
bath and the coupling between qubit and bath. Because 
the coupling between the switched spin is typically of the 
same order as the coupling between the other spins a is 
roughly of the order of 1. 

Second, we consider a solid NMR system, namely 
KPFg. There, we found B m as 90kHz and interspin cou- 
plings ranging from 3.3kHz to 11kHz 29] . This implies 
J/B m as 0.04 — 0.12 whereas a ranges between 0.3 and 
3. Another system is adamantane, for which we assume 
B m as 150kHz and J as 15kHz so that J/B m as 0.1. The 
ratio a is again taken to be of the order of 1 [30] . 

Third, we consider the electronic spin in a quantum 
dot as the qubit. The experimental investigation of tem- 
perature dependent spin relaxation has just started [3l| . 
The pulses are very short (r p as lps) which implies for 
a 7T pulse according to B m T p /h = 0/2 = ir/2 the ampli- 
tude B m as lmeV. But it is much less clear which A or 
a one should consider. In Ref. [3l| a bosonic bath with 
spectral density J c g(uj) is considered. Taking the Debye 
frequency wd = 27.5meV as upper cutoff and deducing 
J from 

J 2 = J e s(u)du (19) 
Jo 

one obtains J ~ 0.3 — 20eV which implies enormous val- 
ues for J/B m but small values for a — uj^ /J. 

But closer inspection of the estimates for T 2 [3l[ re- 
veals that the above estimate is not the relevant one. 
Rather the internal energy scale appears to be set by the 
energy splitting A as 70/ieV of the two qubit states. The 
characteristic coupling is found by restricting the inte- 
gral in ([Tn]) to the interval [0, A]. Then J as 1 — 6neV 
ensues which implies J/B m « 10~ 6 — 10~ 5 and a as 10 4 . 
Hernandez et al. [3l| doubt the relevance of the spin re- 
laxation via Rashba and Dresselhaus terms advocating 
phonon-induced dephasing [3^, [Hj]. Then one should 
rather estimate J 2 as TA with T as 0.2/ieV implying 
J as 4/ieV. Then J/B m as 0.004 and a as 20. This ex- 
ample illustrates that the unambiguous identification of 
the relevant processes of decoherence is still a challenging 
task. 

Fourth, we consider a qubit realized by charge states 
in a superconducting device [341 ] . The pulse length is 
T p as 80ps implying B m as 15/ieV. The coupling J is 
taken from the free decay J as /i/150ps as 5/ieV while we 
deduce ojf, as 0.2/ieV from the decay of the signal with an 
echo pulse. So J/B m as 0.3 and a as 0.04. 

Fifth and last, we look at trapped ions [35[ for which 
we found pulse lengths in the range of microseconds im- 
plying B m as 1MHz. The coupling to optical modes takes 
values J as 20 - 200kHz so that J/B m as 0.02 - 0.2. Less 
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Figure 5: (Color online) Deviation d as a function of J/B m 
(B m is the maximum amplitude of the pulse) for various 
lengths N of the spin chain at a = 5.0. The data refers 
to the ASY1 pulse. 
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obvious is the relevant internal energy scale a>&. The en- 
ergies of the optical modes in the cavities are fairly high 
between 1 and 40Ghz so that a would range in the order 
of 10 6 . Thus the question arises whether this is really the 
relevant scale or whether the very fast modes average out 
so that a much lower effective scale comes into play. 

The above numbers provide a rough guideline in which 
range today's experiments are done. Surely, more elab- 
orate investigations of the relevant decoherence mecha- 
nisms are called for. 



IV. THE NUMERICAL ANALYSIS 

Remarks on the program The numerical data was ob- 
tained using C++ routines. Many of the matrix calcu- 
lations were realized with the help of the MATPACK- 
package [36[. The exponentials of the matrices were cal- 
culated using routines adapted from EXPOKIT [37[ ab- 
breviated padm. These are techniques based on Pade 
summation. Note that this approach is well-suited to 
deal with piecewise constant pulses whereas continuously 
varying pulses are not accessible. 

As anticipated from the analytical calculation, only 
minor finite-size effects occur. This is illustrated numer- 
ically in Fig. \5\ for one particular pulse. But all other 
pulses show the same behavior. Indeed, the finite-size 
effects are completely negligible in the region of small 
values of J. Hence we conclude that a moderate number 
of bath spins is sufficient. In the data presented here we 
routinely use N = 7 and N — 10. For these system sizes 
no particular matrix algorithms are needed. 

7T pulses: vanishing linear order We consider sym- 
metric and asymmetric pulses with angle 9 = 7r which 
satisfy rjn = 7712 = as defined in ([5]). For comparison, 



Figure 6: (Color online) Case of n pulses. The deviation d is 
plotted as function of J/B m for N = 10 and various values of 
a. For an unbiased comparison of the pulses, J is normalized 
to the maximum amplitude B m of the pulses. The notation 
for the pulses refers to Tab. HI The dashed lines ease the 
comparison with pure power laws. 
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Figure 7: (Color online) Case of it pulses. Plot of the pref- 
actors a a in d — a a J 2 + C(J 3 ) for N — 7. The solid lines 
represent the analytical prediction in Eq. (|18|) . 



also the standard pulse with constant amplitude and fi- 
nite is computed. 

Fig.[S]shows the behavior of the deviations d as a func- 
tion of J / B m for representative values of the parameter 
a. Here B m is the maximum amplitude of the pulse. At 
first thought, a plot as function of Jt p appears reason- 
able. But the comparison as function of J/B m is fairer 
because the simple pulses, for instance the standard one 
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SGLPi, need only a smaller amplitude. Hence they can 
experimentally be realized with a shorter duration r p if 
the apparatus restricts the maximum applicable ampli- 
tude. This advantage is accounted for by the plot versus 
J/B m . 

The quadratic behavior of ASYPi and UPi proves that 
the first order corrections are completely cancelled. This 
is not the case for SGLPi for which the numerical data 
display a linear behavior for small J. For large values of 
a, d starts to deviate from the desired quadratic behavior 
even at relative small values of J. This indicates that the 
internal energy scale ujb = otJ becomes important. 

The comparison between the standard pulse SGLPi 
and the optimized ones ASYPi and UPi shows that a 
crossover takes place. For low values of J the pulses with 
vanishing first order outperform the standard pulse due 
to their steeper decrease. At larger values of J the more 
complicated structure of the optimized pulses does not 
pay anymore and SGLPi is slightly better. Note that the 
value of J where the crossover takes place depends on 
the value of a. For low values of a ASYPi and UPi pay 
up to much larger values of J than for large values of a. 

Data such as presented in Fig. [5] is used to determine 
the prefactors a a defined in 



d = a a J 2 + C(J 3 ) 



(20) 



by fits. The fits are made only within the range of valid- 
ity of the quadratic behavior. The results are plotted in 
Fig. El They agree perfectly with the analytical predic- 
tion from Eq. I|18p. For the quantitative comparison the 
coefficients 7721, ^22? and 7723 are explicitly computed for 
ASYPi and UPi2 by means of Eqs. ||SJ), see also Tab. Q] 

7r/2 pulses: vanishing linear order We consider sym- 
metric and asymmetric pulses with angle 9 = n/2 which 
satisfy r/n = r/12 = as defined in ([8]). For comparison, 
also the standard pulse with constant amplitude and fi- 
nite r/ 1 ) is computed. 

Fig.[5]shows the behavior of the deviations d as a func- 
tion of J/B m for representative values of the parameter 
a. Again, the comparison as function of J / B m is fairer 
for the above mentioned reasons. 

The quadratic behavior of ASYPi2 and UPi2 proves 
that the first order corrections are completely cancelled. 
This is not the case for SGLPi2 for which the numerical 
data displays a linear behavior for small J. For large val- 
ues of a, d starts to deviate from the desired quadratic 
behavior even at relative small values of J. This indi- 
cates that the internal energy scale lui, — aJ becomes 
important. 

The comparison between the standard pulse SGLPi2 
and the optimized ones ASYPi2 and UPi2 shows that a 
crossover takes place. For low values of J the pulses with 
vanishing first order outperform the standard pulse due 
to their steeper decrease. At larger values of J the more 
complicated structure of the optimized pulses does not 
pay anymore and SGLPi2 is slightly better. Note that 
the value of J where the crossover takes place depends 
on the value of a. For low values of a ASYP2i and UPi2 
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Figure 8: (Color online) Case of tt/2 pulses. The deviation d is 
plotted as function of J/B m for iV = 10 and various values of 
a. For an unbiased comparison of the pulses, J is normalized 
to the maximum amplitude B m of the pulses. The notation 
for the pulses refers to Tab. [Il] The dashed lines ease the 
comparison with pure power laws. 
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Figure 9: (Color online) Case of n/2 pulses. Plot of the pref- 
actors a a in d — a a J 2 + 0(J 3 ) for N — 7. The solid lines 



represent the analytical prediction in Eq. (JT8 



pay up to much larger values of J than for large values 
of a. 

Note that the gain of the optimized pulses over the 
standard pulse is most significant for low values of a, i.e., 
for a slow internal bath dynamics, ft is less significant for 
fast internal bath dynamics corresponding to large values 
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Figure 10: (Color online) Case of n/2 pulses with partly van- 
ishing quadratic order. The deviation d is plotted as function 
of J/B m for TV = 10 and various values of a. For an unbiased 
comparison of the pulses, J is normalized to the maximum 
amplitude B m of the pulses. The notation for the pulses refers 
to Tab. UU The dashed lines ease the comparison with pure 
power laws. 



of a. 

Data such as presented in Fig. [8] is used to determine 
the prefactors a a defined in (f2"0"|) by fits. The fits are made 
only within the range of validity of the quadratic behav- 
ior. The results are plotted in Fig. [5J They agree per- 
fectly with the analytical prediction in Eq. (p~8|) . For the 
quantitative comparison the coefficients 7721, V22, and 7723 
are explicitly computed for ASYPi2 and UPi2 by means 
of Eqs. ©, see also Tab. QTJ 

The errors in the fits of the prefactors for the 7r/2 pulses 
were determined by hand for some randomly chosen data 
points using various fitting ranges and fitting functions 
such as b J + a a J 2 or a a J 2 + bJ 3 . This analysis provides 
the error estimates of 12% at a = 0.03 to 32% at a = 28 
for pulse UPi2 and about 10% for all a looking at pulse 
ASYPi2. 

7r/2 pulses: partly vanishing quadratic order In the 
previous work in Refs. (2p| we have proven rigorously that 
no 7T pulse can satisfy the second order Eqs. (|8cl8dl8ep . 
For 7r/2 pulses no such proof is known to us. But we were 
not able to find a solution to all five equations Eqs. (O 
either. 

We managed, however, to find solutions which make 
the first four equations (|8al8bl8cl8dp vanish. The advan- 
tage is that the first order vanishes completely and that 
in second order all the terms of order Ao^t 2 vanish also. 
Only the term of order X 2 t 2 persists. We expect such 



pulses, see Fig. [3] and Tab. QTJ to be advantageous for 
systems where the coupling A between qubit and bath is 
very small, but the internal bath dynamics LOb is not. 

Here we propose two possible examples of 7r/2 pulses, 
symmetric and asymmetric, for which 7711 = 7712 = 0, 
7721 = V22 = 0, but 7723 7^ 0. From the above argu- 
ments, we expect that for large values of a, i.e., fairly 
fast baths, the deviation d(J) displays cubic behavior 
at least in some intermediate range. Fig. [TUJ provides 
the corresponding data. Indeed, one clearly identifies an 
intermediate range where cubic behavior is seen. This 
range is fairly small for small values of a (upper panel 
in Fig. [TO)) but grows upon increasing a (middle panel in 
Fig. [TU)) . For the large values of a analyzed in the lower 
panel in Fig. [TUJ the quadratic behavior below the cubic 
range is not even discernible. But we know from Eq. (|18|) 
that is exists. 

We conclude that even a partial vanishing of the second 
order can be very helpful. This conclusion is supported 
by the comparison to data for ASYPi2 and UPi2 which 
have a vanishing first order, but no vanishing second or- 
der terms. As to be expected, we find that for low val- 
ues of J the pulses S2NDPi2 and A2NDPi2 outperform 
ASYPi2 and UPi2. For larger values of J a crossover 
takes place and there is no need to resort to the more 
complicated pulses S2NDPi2 and A2NDPi2. There, all 
pulses behave very much alike. 

Note that the crossover takes place for lower values of 
J if a is large and viceversa for larger values of J if a is 
small. This is related to the fact that the range of cubic 
behavior occurs at larger values of J for small a. For 
large a the range is larger, but shifted to smaller values 
of J. 



V. CONCLUSIONS 

We numerically simulated the effect of designed short 
control pulses on a qubit coupled to a bath of spins. The 
effect of the short pulse can be approximated in leading 
order of the pulse duration r p as a 8 peak. For finite t p , 
however, corrections occur which we know from previous 
analytical calculations. The aim of the present work was 
two-fold. First, we wanted to confirm the analytical re- 
sults by numerical calculations. Second, we intended to 
analyze to which extent the analytically neglected higher 
orders matter. Put differently, we wanted to see whether 
pulses, which are fine-tuned to make the leading correc- 
tions vanish, outperform the standard pulses. 

The numerical results confirm the analytical results in 
all points. The fine-tuned pulses display qualitatively dif- 
ferent power laws in the deviation d as function of Jt p . 
This deviation measures the difference between the ideal 
pulse, multiplied with the evolution due to free decoher- 
ence, and the realistic pulse. For standard pulses, one 
has d oc J. For the fine-tuned pulses we achieve d oc J 2 . 

In restricted parameter ranges, we obtained even d oc 
J 3 for pulses which make certain parts of the second order 
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corrections vanish. Such pulses were not yet discussed 
before. They are only possible for 9 ^ w. 

The second goal has also been achieved by the analysis 
of the real performance in case of the coupling to a spin 
bath. We could show that the fine-tuned pulses outper- 
form the more standard ones in a large range of parame- 
ters. Furthermore, we estimated the relevant parameters 
for a number of generic experiments. These estimates 
show that many experimental setups are such that the 
fine-tuned pulses should improve on the standard pulses. 
But more investigations, both theoretical and experimen- 
tal, are needed to obtain a complete understanding of the 
important decoherence mechanisms. 

For the above reasons we suggest that the choice of 
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